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Abstract. We propose a scheme to prepare nanomechanical oscillators in 
nonclassical steady states, characterized by a pronounced negative Wigner function. In 
our optomechanical approach, the mechanical oscillator couples to multiple laser driven 
resonances of an optical cavity. By lowering the resonance frequency of the oscillator via 
an inhomogeneous electrostatic field, we significantly enhance its intrinsic geometric 
nonlinearity per phonon. This causes the motional sidebands to split into separate 
spectral lines for each phonon number and transitions between individual phonon Fock 
states can be selectively addressed. We show that this enables the preparation of the 
nanomechanical oscillator in a single phonon Fock state. Our scheme can for example 
be implemented with a carbon nanotube dispersively coupled to the evanescent field 
of a state of the art whispering gallery mode microcavity. 



PACS numbers: 85.85.+j,42.50.Dv,42.50.Wk,03.65.Ta 
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1. Introduction 

The manipulation of mechanical degrees of freedom by light scattering, presently 
known as optomechanics, has attracted increasing interest in recent years [U [2]. 
Current research in this direction is largely driven by substantial progress in cooling 
mechanical oscillators to their ground states [31 IU El El IE] and performing displacement 
measurements close to the standard quantum limit [91 [TOj EEE] • These developments set 
the stage for investigations of the quantum regime of macroscopic mechanical degrees 
of freedom, an endeavor that could lead to experimental tests of potential limitations 
of quantum mechanics [12J. Here, we introduce an optoelectromechanical technique to 
prepare a mechanical oscillator in a nonclassical steady state with a negative Wigner 
function by enhancing the anharmonicity of its motion. 

As has been known since the initial formulation of quantum mechanics, the 
dynamics of a purely harmonic quantum system is hard to distinguish from its classical 
counterpart [T3]. As an example, electrical circuits in the superconducting regime do not 
display quantum behaviour unless they feature a nonlinear element such as a Josephson 
junction. Along these lines nanomechanical oscillators that are coupled to nonlinear 
ancilla systems have been considered [HI [T5J Ej. Here, in contrast, we introduce a 
framework that only makes use of intrinsic properties of the mechanical oscillator. This 
approach allows one to prepare the latter in nonclassical states that are stationary states 
of its dissipative motion. The lifetime of these states is thus not limited by the lifetime 
of the oscillator's excitations. 

In fact, doubly clamped mechanical resonators such as nanobeams naturally feature 
an intrinsic, geometric nonlinear contribution to their elastic energy that gives rise to 
a Duffing nonlinearity [T6l [T71 [T8l [19] and typically only becomes relevant for large 
deflection amplitudes. To investigate quantum effects of mechanical degrees of freedom 
[201 |2~T1 122] it is thus of key importance to devise means for enhancing the nonlinearity 
of nanomechanical structures. 

Here we consider inhomogenous electrostatic fields to decrease the harmonic 
oscillation frequency of a nanomechanical oscillator [23]. This procedure enhances 
the amplitude of its zero point motion and therefore leads to an amplification of its 
nonlinearity per phonon, which eventually becomes comparable to the optical linewidth 
of a high finesse cavity. When combined with an optomechanical coupling to the 
cavity, this approach opens up a variety of possibilities for manipulating the mechanical 
oscillator by driving the cavity modes with lasers. In particular the motional sidebands 
split into separated lines for each phonon number n, of. figure [TJ This feature allows 
to selectively transfer population from phonon Fock state \n) to |n+ 1) (\n') to \n' — 1)) 
by driving the spectral line corresponding to n (n') phonons in the blue (red) sideband. 

To demonstrate the capabilities of our approach, we show that suitably adjusted 
laser fields prepare the mechanical oscillator in a stationary state close to a phonon 
Fock state, i.e. in a state for which the quantum behaviour of an oscillator is most 
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Figure 1. Level diagram of the coupled cavity-resonator-system in a shifted 
representation. |n m lj) denotes a single photon state of cavity mode j with n phonons 
in the mechanical mode, while |0) c denotes the vacuum (i.e. the steady state of the 
decoupled cavity). Due to the nonlinearity, the mechanical energy level spacings 
depend on n and can be addressed individually. Only resonant processes are sketched 
with a blue detuned laser acting on the transition |0) m — >■ |l) m and red detuned lasers 
acting on transitions |2) m — > |l) m and |3) m — > |2) m . 

obvious [24J. The highly nonclassical nature of this state manifests itself in a Wigner 
function that reaches the lowest possible value of —2/ it. This can be verified via the 
power spectrum of an optical probe field [251 126] . Here, the stationary nature of the 
motional states considerably simplifies their measurement as compared to phonon Fock 
states that might be prepared by projective measurements [271 E3 129] . 

Our scheme can, for example, be implemented with state of the art 
nanoelectromechanical systems (NEMS) comprising carbon-based mechanical resonators 
[3T1 130] dispersively coupled to high finesse toroidal [32] or fiber based [33] microcavities 
(see figure [2]). Carbon-based oscillators are particularly well-suited as they combine very 
small transverse dimensions with ultra-low dissipation [31] . 

The remainder of the paper is organized as follows. In section 2, we discuss the 
Hamiltonian of the oscillating nanobeam, including its geometric nonlinearity. In section 
3, we then introduce the optomechanical coupling of the beam to a high finesse cavity 
and present the master equation that describes the dynamics of our system including 
the relevant damping mechanisms. In the following section 4, we provide an analysis of 
how driven cavity modes can drive the nanomechanical resonator into a non-classical 
steady state by adiabatically eliminating the photon degrees of freedom. In the sequel 
we analyse important aspects of an implementation in more detail. In section 5 we 
describe the electrostatic technique that lowers the harmonic oscillation frequency of 
the beam and hence increases its nonlinearity per phonon. The following section 6 
describes a possible feasible experimental setup including a carbon nanotube coupled to 
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a microtoroidal cavity. We then present examples of states that could be prepared in 
this setup in section 7. Finally, we discuss a technique to measure the prepared state 
(section 8) and close the paper with our conclusions and an outlook. 

2. Oscillating nanobeam 

The dynamics of the fundamental flexural mode of a thin beam can be described by the 
Hamiltonian [35] , 

P 2 1 8 
H m = — + -m*ulX 2 + -X 4 , (1) 
m 2m, 2 m 4 K J 

where X is the deflection, P its conjugate momentum, u m the mode frequency and m* 

the effective mass. We assume clamped-clamped boundary conditions and apply thin 

beam theory [36]. This yields Co> mj o = c s k(^P) 2 for the mode frequency without applying 

elastic stress or other external forces, and 8 = 0.060 m^u^ k~ 2 for the nonlinearity 

which results from the streching induced by the deflection. Here, L is the suspended 

length, c s is the phase velocity of compressional phonons and k is given by the relevant 

transverse dimension of the beam times a geometric prefactor. 

The single-mode model ([I]) is valid for sufficiently low phonon occupancies and weak 

nonlinearities such that phonon-phonon interactions involving higher harmonics can be 

neglected. Thus we focus on the latter regime and write H m in terms of phonon creation 

(annihilation) operators (b) to obtain 

H m = hcu m tfb + h^(tf + b) 4 (2) 

with nonlinearity A = I^zpm/^ where xzpm = \/h/2m^uj m is the amplitude of the 
oscillator's zero point motion. 

Note that A oc u~ 2 can be significantly enhanced by "softening" the mode, 
i.e. lowering u m . Whereas this can be achieved by applying compressive stress along the 
beam [351 EZ] , we follow here an alternative approach to allow for better tunability. We 
consider an in-plane flexural resonance (cf. figure |2}d) and use tip electrodes to apply 
a strong inhomogeneous static electric field which polarizes the beam. As a dielectric 
body is attracted towards stronger fields, this adds an inverted parabola to the potential 
that counteracts the harmonic contribution due to elasticity. Whence, as will be borne 
out in detail in section 5, the oscillation frequency is reduced [3B1 [231 ES]- 

3. Optomechanics 

We consider a typical optomechanical setup where the displacement X = xzpm 

{b + &t) 

of the nanomechanical oscillator dispersively couples to several laser driven cavity 
resonances. In a frame rotating with the frequencies of the laser fields u^j the system 
Hamiltonian reads 

+ H m , (3) 



At ^7 



a'j + dj 



+ G f + b) 
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where Aj = ojlj —Uj is the detuning of a given laser j from the closest cavity mode with 
frequency Uj and photon annihilation operator aj, VLj/2 = ■sj Pmjftex/ fvJJ-Lj is the driving 
strength associated to that laser with input power P in j and K ex is the corresponding 
external decay rate. The coupling strength Goj = corresponds to the shift of the 
respective cavity resonance per oscillator's deflection [40J. 

We incorporate the damping of the cavity field at rate k and mechanical dissipation 
at rate 7 m in a master equation and shift the cavity normal coordinates to their steady 
state value aj — > aj + aj, with 

aj ~ 2Aj + m ' ^ 
Neglecting higher order terms in the photon-phonon coupling for the regime where 
{a\aj) <C \otj\ 2 [31111110] leads to the linearized Hamiltonian 



+ H m , (5) 



where g m j = 2oijXzpy\G§j is the enhanced optomechanical coupling. The term 
~ |a| 2 (6^ + 6), causing a static shift of the mechanical mode due to radiation pressure, is 
omitted here since it can be compensated by a suitable choice of electrostatic softening 
fields (see section 5). The master equation for the shifted system now reads 

where 

T> C jp = 2ajpa^ — a^ajp — pa^aj , (7) 
V mP = (n + 1) {2bpb ] - b ] bp - pb ] b) 

+ n{2b ] pb-bb ] p- pbtf) , (8) 

with the Bose number n = 1/ [exp(hu m / k^T) — 1] at environment temperature T. In 
equation ^ we have applied a rotating wave approximation that is justified for the high 
mechanical quality factors we are interested in |%D] . 

Equation ^ describes the dynamics of our system. Yet, the behaviour of the 
mechanical mode becomes more apparent after adiabatically eliminating the photon 
degrees of freedom, as we do in the next section. 



4. Reduced master equation 

For \g m j\ C Kwe now derive a reduced master equation (RME) for the mechanical 
motion by adiabatically eliminating the cavity degrees of freedom [2U]. We start by 
applying a rotating wave approximation (RWA) to the mechanical Hamiltonian which 
yields, 

H' m = hu' m b ] b + h^b ] b ] bb, (9) 
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with the shifted frequency u' m = uj m + A' and the nonlinearity per phonon A' = 6A. 
Here, the RWA is warranted provided the relevant phonon numbers n fulfil n 2 \' <C Qu m . 
Within the RWA, the mechanical energy eigenstates are Fock states and the large 
nonlinearity, A' ^> \g m j\ 2 /K, n7 m , we are interested in naturally leads to Fock number 
resolved dynamics. In an interaction picture with respect to H' m transition operators 
between neighbouring phonon Fock states rotate at frequencies corresponding to the 
respective mechanical energy level spacings, as 



c 



j/H be -iH^t/H = e- iSnt b n , (10) 



with b n = ^Jn\n — l)(n| and 

S n = (E n - E n ^)/h + X'(n - 1) . (11) 

We further focus on the resolved sideband regime for the softened frequency n <C LJ m , 
which is suitable for efficient back action cooling, as well as <C |<7mj| 2 / K to 

ensure that the influence of the cavity modes on the mechanical motion overwhelms 
environmental heating. The resulting RME for the reduced density operator of the 
nanomechanical resonator, /i = Tr c {p}, reads 



l 

A = - 



h 



An 

+ EE -^Vy + ^V^, (12) 



n,j a- 



where the driven cavity modes act like an additional bath for each number transition 
with damping terms 

Z>> = 2byb n -b n bltt-fib n bl, (13) 
V n _ii = 2h ni ibl-b\b n[ i-[ib\b n . (14) 

Note that the frequency shifts 

A (n) = 2 { Aj + Sn + Aj - 6n X (15) 

^ 9 ^ \ 4( A, + 5 n f + k) + 4(A, - 5 n y + K]j 1 ° J 

are negligible compared to u' m , A'. The cavity-induced transition rates between 
neigbouring Fock states 

I 1 2 
A n _ \9m,j\ K 

however, can exceed the mechanical damping A± ■ > / -y m n which leads to a significant 
impact on the mechanical steady state. 



The steady state of the RME ( 12 ) is diagonal in the Fock basis 



fx = J2 P n\n)(n\. (17) 



Hence the relations 

Pn E 7 ^+,i+7mn 



P n _l £ • A n j + 7m [n + 1] 
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determine all occupation probabilities P n and therefore the mechanical state /i. The 



relations (18) show the central feature of the driven nonlinear nanoresonator. Here 
|flW| 2 /« > 7m^ implies P n /P n -i « For a detuning Aj = ±5 n 

we find A± j > A™ • and the corresponding laser drive therefore enhances (decreases) 

Pn/ Pn-l- 

This feature can be used to approach phonon Fock states |n), as illustrated in 
figure [T] for n = 1. Fock states are highly nonclassical and for n odd they exhibit 
a maximally negative Wigner function at the origin. Ensuring that each laser drive 
preferentially addresses only one transition |n) — > \n ± 1) naturally requires A' > k. 
This can be achieved by electrostatic mode softening as we describe in the next section. 

5. Dielectric softening 

The harmonic frequency u m of the nanoresonator can be reduced by adding an 
electrostatic potential V es oc —X 2 to the elastic potential, where X is the deflection 
of the resonator. This is achieved by placing charged electrodes at both sides of the 
nanoresonator (see figure [2]), that generate a strongly inhomogenous electrostatic field. 
To estimate the potential that softens the beam, we consider its electrostatic energy line 
density 

W(x,y) = -I[ a ,|Sf(ar,y) + a ± E 2 ± (x,y)] . (19) 

Here, y and x are the co-ordinates along the beam and the direction of its deflection, 
E\\,i_ are the electrostatic field components parallel and perpendicular to the beam and 
aiu the respective screened polarizabilities. The electrostatic energy as a functional of 
the deflection x(y) thus reads 

V es = [ W{x{y),y)dy. (20) 



We expand W(x(y),y) to second order in the transverse deflection x(y) from its 
equilibrium configuration x(y) = 0. The leading contributions for the fundamental 
mode with profile (f>o(y) can be found by considering x(y) = (fio(y)X. They read 

V es « Ks,o + V es>1 X + V eSt2 X 2 , (21) 

where V es fl is an irrelevant constant and 

Ves, = £ dy d -^^\ x = My) > (22) 

v es , 2 = 2j o d y — — L=oMy)- (23) 

The linear term V es ^X shifts the beam's equilibrium position and can be employed to 
compensate shifts induced by the dispersive coupling to the cavity fields by satisfying 
Ks,i = Go\®j\ 2 - The quadratic term V CSt 2X 2 in turn is negative and reduces the 

fundamental frequency as 

u4 = u£ j0 --^|V es , 2 |. (24) 
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Figure 2. a) Possible setup with a carbon nanotube dispersively coupled to the 
evanescent field of a microtoroid and tip electrodes used for the electrostatic softening 
of the mechanical frequency (blue cylinders), b) NEMS chip carrying the nanotube 
and the electrodes. The oscillation direction of the nanotube is indicated by the red 
arrows. In a), nanotube and electrodes are not drawn to scale with the toroid and the 
NEMS chip is not shown for illustration purposes. The placement of nanotube and 
electrodes minimizes degradation of the cavity finesse and maximizes optomechanical 
coupling. 

At | V^ S; 2 1 = o the buckling instability [35], [37] occurs and by tuning the electrostatic 

fields below it one can reduce u m by any desired factor ( = w mi0 /w m . Note that the 
softening increases the zero point motion amplitude and, hence, also the optomechanical 
coupling per photon Go^zpm- A particularly useful feature of this dielectric approach 
to soften the nanobeam is that the electrostatic potential and hence u m , Go^zpm and 
A can be tuned in situ. To illustrate the feasibility of our approach we now turn the 
description of a feasible experimental setup for its implementation. 



6. Carbon nanotube and microtoroid 



As a suitable device for realizing our scheme we envisage a setup with a single-walled 
carbon nanotube as the mechanical oscillator (c s = 21000ms -1 ) that interacts with the 
evanescent field of a whispering gallery mode of a microtoroid cavity (32], cf. figure. [2] 
A carbon nanotube is a favorable candidate since without softening (u m = u) mj o) one 
obtains A oc (m*R 2 )~ l indicating that low effective mass m* and small transverse 
dimension oc h are desirable. For a carbon nanotube of radius R, one has R = R/V2. 
We consider a wide band gap semiconducting carbon nanotube which for the infrared 
laser- wavelengths involved will behave as a dielectric [UJ. Taking into account that the 
transverse optical polarizability is negligible (i.e. depolarization) and approximating the 
longitudinal one by its static value, the optomechanical coupling can be estimated as 

(cf. ro 

G o,j ~ ujj^^K^-^C, (25) 
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where V c = 7ia^L c is the mode volume of the whispering gallery mode with cavity 
length L c and effective mode waist 2a c . £ is the ratio of the field at the cavity 
surface to the maximum field inside, d the distance between the chip carrying the 
nanotube and the microtoroid's surface, and k'J 1 the decay length of the evanescent 
field. The correction factor C depends on the type of whispering gallery mode and 
the placement of the carbon nanotube. We assume for simplicity that the NEMS and 
photonic chips are perpendicular to each other and find for a TE mode and optimal 
placement C ~ 0.17/ ^/ K±(d + a c ). This optimum involves a small displacement of the 
carbon nanotube midpoint from the equatorial plane of the toroid to ensure coupling to 
the in-plane motion and an angle of ~ 50° between the carbon nanotube axis and the 
equatorial plane of the toroid to circumvent depolarization. 

In this setup, the tip electrodes used to lower the mechanical oscillation frequency 
need to be positioned close to the cavity surface to ensure a sufficiently inhomogenous 
electrostatic field while maintaining the optomechanical coupling high enough. Hence, 
a concomitant degradation of the cavity's finesse due to photon scattering or absorption 
at the tip electrodes is a natural concern. We find that cavity losses are minimized if one 
uses TE whispering gallery modes and tip electrodes of sub-wavelength diameter d e \ that 
are aligned approximately parallel to the equatorial plane of the toroid (cf. figure. [2^i 
and|2}3). For this arrangement the dominant loss mechanism is absorption. The ratio 
of absorbed power P a to power circulating in the cavity P c can be estimated to be (cf. 
[42] . chapter 8), 



where a is the 2D optical conductivity of the electrodes and 6 their misalignment with 
respect to the equatorial plane of the toroid. The resulting finesse is given by 



where J- c is the bare cavity finesse. For relevant cavity parameters we find <C J 7 ^ 1 
if we assume 6 < 1°, d e \ < 10 nm and a < 2 x 10~ 5 Q . Metallic carbon nanotubes, for 
example, fulfil the latter two requirements provided their resonances do not match the 
relevant cavity frequencies Uj [43J. 

In the following we now give some examples of steady states for the mechanical 
oscillator that can be generated in the setup described here. 

7. Examples 

To illustrate the potential of our scheme, we show in figure [3] results for a (10, 0) carbon 
nanotube characterised by R = 0.39 nm, ot\\ = 142(47re A 2 ) and a± = 10.9(47re A 2 ) 04]. 
The nanotube has a length of L = 1.0 /im and is coupled to whispering gallery modes 
with wavelengths 2^c/ojj ~ 1.1 /im, a c = 1.4/im, and £ = 0.4 of a silica microtoroid 
with finesse J- c = 3 x 10 6 [45J and circumference L c = 1.35 mm. The applied voltages 




(26) 




(27) 
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Figure 3. Steady state of an oscillating (10,0) carbon nanotube with L = 1.0 (im 
and Q m = 5 x 10 6 for T = 20 mK, J c = 3x 10 6 and L c = 1.35 mm. a) Fock state 
occupation probabilities P n . b) Wigner function W . These results are obtained from 
a numerical solution of equation Q6u . 



are chosen so that tip electrodes at a distance of 20 nm from the tube induce maximal 
fields at its axis of E\\ « 1.20 x 10 7 V/m and E± w 1.778 x 10 6 V/m and lower its 
mechanical frequency by a factor £ = 4.0. For the optomechanical couplings, we consider 
g m> j/27T = 21.0 kHz for all j which can be achieved with d = 50 nm, K ex /n = 0.1 and 
a launched power P inj = 1.2 W. The above parameters lead to u) m /2ir = 5.23 MHz, 
\/2n = 209 kHz and k/2tt = 52.3 kHz. We assume carbon nanotube Q-values of 
Qm = 5 x 10 6 and Q m = 5 x 10 5 [fjand an environmental temperature T = 20 mK. 

For these parameters and for laser detunings Aj ~ 8\ and A 2 3 —£2,3 a single 
phonon Fock state can be prepared with 91% fidelity (for Q m = 5 x 10 6 ) as shown in 
figure [3] This steady state of the mechanical oscillator has a Wigner function that is 
clearly negative at the origin, W(0,0) = -0.53 for Q m = 5 x 10 6 and W(0,0) = -0.30 
for Q m = 5 x 10 5 . 

The results shown in figure [3] have been obtained from a numerical solution of 
equation ^ with numerically optimised detunings Aj. These numerical results do not 
rely on the rotating wave approximation used in equation rt9J). Here the fastly rotating 
terms lead to corrections to the energy level spacings 5 n that can be comparable to k. 
Hence to match the laser detunings with sufficient accuracy to the energy level spacings 
8 n , the latter have been determined numerically. 

In turn, for the more moderate parameters: L = 1.7 fim, T c = 2 x 10 6 , L c = 1.80 mm 
(k/2ix = 64.2kHz), ( = 3.3, (which imply oo m /2n = 2.13MHz, and \/2ix = 85.5kHz), 
launched laser powers Pm 5 i/2 = 22 mW and Pi ni 3 = 44 mW, Q m = 1.5 x 10 6 , and 
T = 30 mK, we still find a significant negative peak of depth W min = —0.15 (see figure [4]). 

I For carbon nanotubes ^-Q m ~ 10 14 Hz has already been achieved (G.A. Steele, private 
communication) and there are strong indications that mechanical dissipation in doubly clamped CNTs 
strongly depends on the amplitude of the motion |49) . Hence even larger Q- values would result for the 
amplitudes comparable to the zero point motion that are relevant here. 




Figure 4. Steady state of an oscillating (10,0) carbon nanotube with L = 1.7 /mi 
and Q m = 1.5 x 10 6 for T = 30 mK, T c = 2 x 10 6 and L c = 1.8 mm. a) Fock state 
occupation probabilities P n . b) Wigner function W . These results are obtained from 
a numerical solution of equation Q6u . 



We note that the deep resolved-sideband and undercoupled conditions used here, 
where the frequencies of the driving lasers are detuned by several cavity linewidths 
from the optical resonances of the cavity and losses into the coupling fiber are only a 
small fraction of the total losses, imply that the power absorbed inside the cavity is 
much smaller than the launched power (below 30 /xW in the example of figure [3] and 
below 7 fiW in the example of figure [1]). Whereas an intracavity absorption ~ 1 /iW has 
already been shown to be compatible with a cryogenic environment |46j . even higher 
finesses and lower absorption could be attained by using crystalline resonators [47J. 
Additionally, the CNT polarizabilities could be enhanced well above the static values 
we have assumed by profiting from the corresponding excitonic optical resonances [4~8] 
allowing for a substantial reduction in the laser powers required. 

Furthermore, in both examples, the voltages applied to soften the frequency of the 
CNT nanomechanical resonator generate electric fields at the CNT that are well below 
its threshold for dielectric breakdown and are comparable to fields applied in recent 
experiments [23] . 

Finally, we note that our approach can be generalized to the preparation of Fock 
states with n > 1. 

8. Measurements 

To verify whether a nonclassical steady state with a negative Wigner function W has 
been successfully prepared, one can measure sidebands of the output power spectrum 
S(oj) of an additional cavity mode weakly driven by a probe laser with frequency Co>l- 
Within our quantum noise approach [ID] and under the same conditions for which Eq. 
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Figure 5. Power spectrum S(lu) of a resonant probe laser with the sidebands split 
into well resolved lines corresponding to the phonon numbers n. This example shows 
the result of equation ( 28 1 for the parameters also used in figure [31 



(12) holds we find for the sidebands (i.e. where u 7^ u£) 



S(u) oc 



nT n A[ 



P-n 



(u — ojl — 5 n 



r 2 
4 



+ 



nV n A^ 



P, 



n-l 



(2f 



with sideband linewidths 

T n = n[Al + A n + + lm (2n + l)} 

+ (n-l)[Ar 1 + 7m(^+l)] 
+ (n + l)[^ +1 + 7m n]. 



(29) 



Since the negativity of the Wigner function is in our scheme expected to be maximal 
at the origin where W(0, 0) = - Yln^~^) n the nonclassicality of the state can be 
read off from S{uS). For a probe laser that is resonant with a cavity mode, one has 
A\ = A™ and the ratios of peak intensities S(ul + S n )/S(ui, — S n ) ~ P n /P„_i together 
with 1 = J2 n P n allow to determine the P n using that occupations above a certain Fock 
number are negligible. 



9. Conclusions and Outlook 



We have shown that a dielectrically enhanced nonlinearity in nanomechanical oscillators 
enables controlling their motion at a Fock state resolved level. The scheme set forth 
opens new avenues for optomechanical manipulations of nanobeams, in particular by 
exploiting the in situ tunability of the electrostatic softening approach introduced. 
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